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Abstract 

We introduce and analyze the physics of "driving reversal" experiments. These are pro- 
totype wavepacket dynamics scenarios probing quantum irreversibility. Unlike the mostly 
hypothetical "time reversal" concept, a "driving reversal" scenario can be realized in a lab- 
oratory experiment, and is relevant to the theory of quantum dissipation. We study both 
the energy spreading and the survival probability in such experiments. We also introduce 
and study the "compensation time" (time of maximum return) in such a scenario. Exten- 
sive effort is devoted to figuring out the capability of either Linear Response Theory (LRT) 
or Random Matrix Theory (RMT) in order to describe specific features of the time evolu- 
tion. We explain that RMT modeling leads to a strong non-perturbative response effect that 
differs from the semiclassical behavior. 
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1 Introduction 



In recent years there has been an increasing interest in understanding the theory 
of driven quantized chaotic systems [1,2,3,4,5,6,7,8,9,10,11]. Driven systems are 
described by a Hamiltonian 7i((5, -P, x{t)), where x{t) is a time dependent param- 
eter and (Q, P) are some generalized actions. Due to the time dependence of x{t), 
the energy of the system is not a constant of motion. Rather the system makes 
"transitions" between energy levels, and therefore absorbs energy. This irreversible 
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loss of energy is known as dissipation. To have a clear understanding of quantum 
dissipation we need a theory for the time evolution of the energy distribution. 

Unfortunately, our understanding on quantum dynamics of chaotic systems is still 
quite limited. The majority of the existing quantum chaos literature concentrates on 
understanding the properties of eigenfunctions and eigenvalues. One of the main 
outcomes of these studies is the conjecture that Random Matrix Theory (RMT) 
modeling, initiated half a century ago by Wigner [12,13], can capture the universal 
aspects of quantum chaotic systems [14,15]. Due to its large success RMT has 
become a major theoretical tool in quantum chaos studies [14,15], and it has found 
applications in both nuclear and mesoscopic physics (for a recent review see [16]). 
However, its applicability to quantum dynamics was left unexplored [17,18]. 

This paper extends our previous reports [10,17,18] on quantum dynamics, both in 
detail and depth. Specifically, we analyze two dynamical schemes: The first is the 
so-called wavepacket dynamics associated with a rectangular pulse of strength +e 
which is turned on for a specified duration; The second involves an additional pulse 
followed by the first one which has a strength — e and is of equal duration. We define 
this latter scheme as driving reversal scenario. We illuminate the direct relevance of 
our study with the studies of quantum irreversibility of energy spreading [10] and 
consequently with quantum dissipation. We investigate the conditions under which 
maximum compensation is succeeded and define the notion of compensation (echo) 
time. To this end we rely both on numerical calculations performed for a chaotic 
system and on analytical considerations based on Linear Response Theory (LRT). 
The latter constitutes the leading theoretical framework for the analysis of driven 
systems and our study aims to clarify the limitations of LRT due to chaos. Our re- 
sults are always compared with the outcomes of RMT modeling. We find that the 
RMT approach fails in general, to give the correct picture of wave-evolution. RMT 
can be trusted only to the extend that it gives trivial results that are implied by per- 
turbation theory. Non-perturbative effects are sensitive to the underlying classical 
dynamics, and therefore the — > behavior for effective RMT models is strikingly 
different from the correct semiclassical limit. 

The structure of this paper is as follows: In the next section we discuss the no- 
tion of irreversibility which is related to driving reversal schemes and distinguish 
it from micro-reversibility which is associated with time reversal experiments. In 
Section 3 we discuss the driving schemes that we are using and we introduce the 
various observables that we will study in the rest of the paper. In Section 4, the 
model systems are introduced and an analysis of the statistical properties of the 
eigenvalues and the Hamiltonian matrix is presented. The Random Matrix Theory 
modeling is presented in Subsection 4.4. In Section 5 we introduce the concept of 
parametric regimes and exhibit its applicability in the analysis of parametric evo- 
lution of eigenstates [19]. Section 6 extends the notion of regimes in dynamics and 
presents the results of Linear Response Theory for the variance and the survival 
probability. The Linear Response Theory (LRT) for the variance is analyzed in de- 
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tails in the following Subsection 6.1. In this subsection we also introduce the notion 
of restricted quantum-classical correspondence (QCC) and show that, as far as the 
second moment of the evolving wavepacket is concerned, both classical and quan- 
tum mechanical LRT coincides. In 6.5 we present in detail the results of LRT for 
the survival probability for the two driving schemes that we analyze. The following 
Sections 7 and 8 contain the results of our numerical analysis together with a crit- 
ical comparison with the theoretical predictions obtained via LRT. Specifically in 
Sect.7, we present an analysis of wavepacket dynamics [18] and expose the weak- 
ness of RMT strategy to describe wavepacket dynamics. In Sect.8 we study the 
evolution in the second half of the driving period and analyze the Quantum Irre- 
versibility in energy spreading, where strong non-perturbative features are found 
for RMT models [10]. Section 9 summarizes our findings. 



2 Reversibility 

The dynamics of either a classical or a quantum mechanical system is generated 
by a Hamiltonian Ti.{Q, P; x{t)) where x = (Xi, X2, X3, ...) is a set of parameters 
that can be controlled from the "outside". In principle x stands for the infinite set 
of parameters that describe the electric and magnetic fields acting on the system. 
But in practice the experimentalist can control only few parameters. A prototype 
example is a gas of particles inside a container with a piston. Then Xi may be the 
position of the piston, X2 may be some imposed electric field, and X3 may be some 
imposed magnetic field. Another example is electrons in a quantum dot where some 
of the parameters X represent gate voltages. 

What do we mean by reversibility? Let us assume that the system evolves for some 
time. The evolution is described by 



where Exp stands for time ordered exponentiation. In the case of the archetype ex- 
ample of a container with gas particles, we assume that there is a piston (position 
X) that is translated outwards increasing). Then we would "undo" the evo- 

lution, by displacing the piston "inwards" (X^{t) decreasing). In such a case the 
complete evolution is described by U[x] = U[x^]U[x"^]. If we getU = 1 (up to a 
phase factor), then it means that it is possible to bring the system back to its original 
state. In this case we say that the process U[x] is reversible. 

In the strict adiabatic limit the above described process is indeed reversible. What 
about the non-adiabatic case? In order to have a well posed question we would like 
to distinguish below between "time reversal" and "driving reversal" schemes. 




(1) 
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2.1 Time reversal scheme 



Obviously we are allowed to invent very complicated schemes in order to "undo" 
the evolution. The ultimate scheme (in the case of the above example) involves 
reversal of the velocities. Assume that this operation is represented by Ut, then the 
reverse evolution is described by 

[/reverse = UtU[x^]Ut , (2) 

where in xsit) we have the time reversed piston displacement (X{t)) together 
with the sign of the magnetic field (if it exists) should be inverted. The question is 
whether Ut can be realized. If we postulate that any unitary or anti-unitary transfor- 
mation can be realized, then it follows trivially that any unitary evolution is "micro- 
reversible". But when we talk about reversibility (rather than micro-reversibility) 
we allow control over a restricted set of parameters (fields). Then the question is 
whether we can find a driving scheme, named x^, such that 

Ut = U[x'^] ??? (3) 

With such restriction it is clear that in general the evolution is not reversible. 

Recently it has been demonstrated in an actual experiment that the evolution of spin 

system (cluster with many interacting spins) can be reversed. Namely, the complete 
evolution was described by U[x] = U[x'^]U[x^]U[x'^]U[x^], where U[x^] is gener- 
ated by some Hamiltonian Ha = Hq + eW. The term Hq represents the interaction 
between the spins, while the term W represents some extra interactions. The uni- 
tary operation U[x'^] is realized using NMR techniques, and its effect is to invert 
the signs of all the couplings. Namely U[x'^]HoU[x'^] — —Hq. Hence the reversed 
evolution is described by 

[/reverse = exp (^-^^(-Tio + £>V)) , (4) 

which is the so-called Loschmidt Echo scenario. In principle we would like to have 
£ so as to get [/ — 1, but in practice we have some un-controUed residual fields 
that influence the system, and therefore e ^ 0. There is a huge amount of literature 
that discusses what happens in such scenario [20,21,22,23,24]. 



2.2 Driving reversal scheme 

The above described experiment is in fact exceptional. In most cases it is possible 
to invert the sign of only one part of the Hamiltonian, which is associated with the 
driving field. Namely, if for instance f/[x^] is generated by Ha — Hq-\- eW, then 
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we can realize 

tA-everse = exp (^--^^(7^0 - , (5) 

whereas Eq.(4) cannot be realized in general. We call such a typical scenario "driv- 
ing reversal" in order to distinguish it from "time reversal" (Loschmidt Echo) sce- 
nario. 

The study of "driving reversal" is quite different from the study of "Loschmidt 
Echo". A simple minded point of view is that the two problems are formally equiva- 
lent because we simply permute the roles of Ho and W. In fact there is no symmetry 
here. The main part of the Hamiltonian has in general an unbounded spectrum with 
well defined density of states, while the perturbation W is assumed to be bounded. 
This difference completely changes the "physics" of dynamics. 

To conclude the above discussion we would like to emphasize that micro-reversibility 
is related to "time reversal" experiment which in general cannot be realized, while 
the issue of reversibility is related to "driving reversal", which in principle can be 
realized. Our distinction reflects the simple observation that not any unitary or anti- 
unitary operation can be realized. 



3 Object of the Study 



In this paper we consider the issue of irreversibility for quantized chaotic systems. 
We assume for simplicity one parameter driving. We further assume that the varia- 
tion of x{t) is small in the corresponding classical system so that the analysis can 
be carried out with a linearized Hamiltonian. Namely, 

niQ,P;x{t))^no + 5x{t)W , (6) 

where Ho = Ti.{Q, P; x{0)) and 6x = x{t) — x{Q). For latter purposes it is conve- 
nient to write the perturbation as 

5xit) = e X fit) , (7) 

where e controls the "strength of the perturbation", while f{t) is the scaled time 
dependence. Note that if f{t) is a step function, then e is the "size" of the pertur- 
bation, while if f{t) (xt then e is the "rate" of the driving. In the representation of 
Ho we can write 

H^E + Sx{t)B , (8) 

where by convention the diagonal terms of B are absorbed into the diagonal matrix 
E. From general considerations that we explain later it follows that S is a banded 
matrix that looks random. This motivates the study of an Effective Banded Ran- 
dom Matrix (EBRM) model, as well as its simplified version which is the standard 
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Wigner Banded Random Matrix (WBRM) model. (See detailed definitions in the 
following). 

In order to study the irreversibility for a given driving scenario, we have to introduce 
measures that quantify the departure from the initial state. We define a set of such 
measures in the following subsections. 

3. 1 The evolving distribution (n | no ) 

Given the Hamiltonian P; x), an initial preparation at state |no), and a driving 
scenario x{t), it is most natural to analyze the evolution of the probability distribu- 
tion 

Pt(n|no) = |(n|C/(t)|no)r (9) 
We always assume that x{t) — a;(0). 

By convention we order the states by their energy. Hence we can regard Pt(n|no) 
as a function ofr — n — riQ, and average over the initial preparation, so as to get a 
smooth distribution Vt{r). 

The survival probability is defined as 

V{t) = |(no|C/W|no)|' = PtKl^o) , (10) 

and the energy spreading is defined as 

5E{t) = lY.PMno){E^-ErJ^ . (11) 
y n 

These are the major measures for the characterization of the distribution. In later 
sections we would like to analyze their time evolution. 

The physics of 5E{t) is very different from the physics of V{t) because the former 
is very sensitive to the tails of the distribution. Yet, the actual "width" of the dis- 
tribution is not captured by any of these measures. A proper measure for the width 
can be defined as follows: 

SEcore{t) = K5% - n25%]A , (12) 

where A is the mean level spacing and Uq is determined through the equation 
X^„Pt(n|no) = q. Namely it is the width of the main body of the distribution. 
Still another characteristic of the distribution is the participation ratio Snip^{t). It 
gives the number of levels that are occupied at time t by the distribution. The ratio 
^^ipr/ (^75% ~ ^25%) can be used as a measure for sparsity. We assume in this paper 
strongly chaotic systems, so sparsity is not an issue and SrijpR ~ 5£'core/A. 
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Fig. 1. Shape of the appHed driving schemes f{t); wavepacket dynamics (left panel) and 
driving reversal scenario (right panel) 



3.2 The compensation time tr 



In this paper we consider two types of driving schemes. Both driving schemes are 
presented schematically in Figure 1. 

The first type of scheme is the wavepacket dynamics scenario for which the driving 
is turned-on at time t — and turned-off at a later time t — T. 

The second type of scenario that we investigate is what we call driving reversal. In 
this scenario the initial rectangular pulse is followed by a compensating pulse of 
equal duration. The total period of the cycle is T. 

In Figure 9 we show representative results for the time evolution of 6E{t) in a 
wavepacket scenario, while in Figure 12 we show what happens in case of a driving 
reversal scenario. Corresponding plots for V{t) are presented in Figure 13. We shall 
define the models and we shall discuss the details of these figures later on. At this 
stage we would like to motivate by inspection of these figures the definition of 
"compensation time". 

We define the compensation time tr, as the time after the driving reversal, when 
maximum compensation (maximum return) is observed. If it is determined by the 
maximum of the survival probability kernel V{t), then we denote it as t^. If it is 
determined by the minimum of the energy spreading SE{t) then we denote it as tf. 
It should be remembered that the theory of V{t) and SE{t) is not the same, hence 
the distinction in the notation. The time of maximum compensation is in general 
not tr — T but rather 

T/2<tr<T . (13) 

We emphasize this point because the notion of "echo", as used in the literature, 
seems to reflect a false assertion [24]. 

For the convenience of the reader we concentrate in the following table on the ma- 
jor notations in this paper: 
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Notation 


explanation 


reference 


n{Q,P;x{t)) 


classical linearized Hamiltonian 


Eq.(6) 


m 


generalized force 


Eq.(17) 


C{t) 


correlation function 


Eq.(18) 




correlation time 


- 


C{u) 


fluctuation spectrum 


Eq.(19) 


n^E + 5xB 


The Hamiltonian matrix 


Eq.(8) 


2DW 


the physical model system 


Eq.(14) 


EBRM 


the corresponding RMT model 


- 


WBRM 


the Wigner RMT model 


- 


A 


mean level spacing 


Eq.(23) 


A, 


energy bandwidth 


Eq.(24) 


a 


RMS of near diagonal couplings 


- 


-fspacings (s) 


energy spacing distribution 


Eq.(16) 


-^couplings (?) 


distribution of couplings 


Eq.(25) 


En{x) 


eigen-energies of the Hamiltonian 


- 


En-E^ ^ rA 


estimated energy difference for r = n — m 


- 


P{n\m) 


overlaps of eigenstates given a constant perturbation e 


Eq.(27) 


P{r) 


smoothed version of P{n\m) 


- 


V{5x) 


the number of levels that are mixed non-perturbatively 


- 


5Ec\ oc 5x 


the classical width of the LDoS 


Eq.(29) 


6x = efit) 


driving scheme 


Eq.(7) 


T 


The period of the driving cycle (if applicable) 


- 


Pt{n\m) 


the transition probability 


Eq.(9) 


Pt{r) 


smoothed version of Pt(n no) 


- 


V{t) 


the survival probability Pt(no no) 


Eq.(lO) 


p{t) = 1 - V{t) 


total transition probability 


Eq.(47) 


5E{t) 


energy spreading 


Eq.(ll) 


^Ecoreit) 


the "core" width of the distribution 


Eq.(12) 




compensation time for the survival probability 


- 




compensation time for the energy spreading 




tpn^ ^sdn) ^erg 


various time scales in the dynamics 


Eq.(59,62,70) 


Spit 


borders between regimes 


Eq.(31,33) 


-fpOPT, -fprt) -Psc 


various approximations to P() 


Eq.(30,32,34,35) 
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Fig. 2. Equipotential contours (left) of the model Hamiltonian Ho for different energies and 
the Poincare section (right) of a selected trajectory at E = 3. Some tiny quasi-integrable 
islands are avoided (mainly at (0, 0)). 

4 Modeling 

We are interested in quantized chaotic systems that have few degrees of freedom. 
The dynamical system used in our studies is the Pullen-Edmonds model [25,26]. It 
consists of two harmonic oscillators that are nonlinearly coupled. The correspond- 
ing Hamiltonian is 

n{Q, P;x) = ^ (P^ + P| + Ql + Ql) + xQlQl . (14) 

The mass and the frequency of the harmonic oscillators are set to one. Without loss 
of generality we set x(0) = xq = 1. Later we shall consider classically small defor- 
mations (5x <^ 1) of the potential. One can regard this model (14) as a description 
of a particle moving in a two dimensional well (2DW). The energy E is the only di- 
mensionless parameter of the classical motion. For high energies E > 5 the motion 
of the Pullen-Edmonds model is ergodic. Specifically it was found that the measure 
of the chaotic component on the Poincare section deviates from unity by no more 
than 10-3[25,26]. 

In Figure 2 we display the equipotential contours of the model Hamiltonian (14) 
with xq — 1. We observe that the equipotential surfaces are circles but as the en- 
ergy is increased they become more and more deformed leading to chaotic motion. 
Our analysis is focused on an energy window around ~ 3 where the motion is 
mainly chaotic. This is illustrated in the right panel of Figure 2 where we report the 
Poincare section (of the phase space) of a selected trajectory, obtained from Hq at 
E — 3. The ergodicity of the motion is illustrated by the Poincare section, filling 
the plane except from some tiny quasi- integrable islands. 

The perturbation is described by W = QlQl- In the classical analysis there is only 
one significant regime for the strength of the perturbation. Namely, the perturbation 
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is considered to be classically small if 

5x < Ea , (15) 

where Scx — l. This is the regime where (classical) linear analysis applies. Namely, 
within this regime the deformation of the energy surface Tio = E can be described 
as a linear process (see Eq. (29)). 

4.1 Energy levels 

Let us now quantize the system. For obvious reasons we are considering a de- 
symmetrized 1/8 well with Dirichlet boundary conditions on the lines Qi = 0, 
(^2 = and Qi = Q2. The matrix representation of Ho in the basis of the un- 
coupled system is very simple. The eigenstates of the Hamiltonian Ho are then 
obtained numerically. 

As mentioned above, we consider the experiments to take place in an energy win- 
dow 2.8 < -E < 3.1 which is classically small and where the motion is predomi- 
nantly chaotic. Nevertheless, quantum mechanically, this energy window is large, 
i.e., many levels are found therein. The local mean level spacing A(£') at this en- 
ergy range is given approximately by A ~ A.3h^. The smallest h that we can 
handle is h = 0.012 resulting in a matrix size of about 4000 x 4000. Unless stated 
otherwise, all the numerical data presented below correspond to a quantization with 
fi = 0.012. 

As it was previously mentioned in the introduction, the main focus of quantum 
chaos studies has so far been on issues of spectral statistics [14,15]. In this context 
it turns out that the sub -h statistical features of the energy spectrum are "universal", 
and obey the predictions of RMT. In particular we expect that the level spacing 
distribution P(s) of the "unfolded" (with respect to A) level spacings s„ = {En+i — 
En)/ A will follow with high accuracy the so-called Wigner surmise. For systems 
with time reversal symmetry it takes the form [14,27] 

^spacings (s) = I S e" ? (16) 

indicating that there is a linear repulsion between nearby levels. Non-universal (i.e. 
system specific) features are reflected only in the large scale properties of the spec- 
trum and constitute the fingerprints of the underlying classical chaotic dynamics. 

The de-symmetrized 2DW model shows time reversal symmetry, and therefore we 

expect the distribution to follow Eq.(16). The analysis is carried out only for the 
levels contained in the chosen energy window around E = 3. Instead of plotting 
P(s) we show the integrated distribution I{s) = Jq P{s')ds', which is independent 
of the bin size of the histogram. In Figure 3 we present our numerical data for /(s) 
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Fig. 3. The integrated level spacing distribution I{S) of the unperturbed Hamiltonian Tio 
(h = 0.012). The dashed line is the theoretical prediction for the GOE. Inset: Difference 
between the theoretical prediction /th(>S') and the actual distribution I{S). 

while the inset shows the deviations from the theoretical prediction (16). The agree- 
ment with the theory is fairly good and the level repulsion is clearly observed. The 
observed deviations have to be related on the one hand to the tiny quasi-integrable 
islands that exist at = 3 as well as to rather limited level statistics. 



4.2 The band-profile 

In this subsection we explain that the band-structure of B is related to the fluctua- 
tions of the classical motion. This is the major step towards RMT modeling. 

Consider a given ergodic trajectory {Q{t),P{t)) on the energy surface 
'H{Q{0), -P(O); xo) = E. An example is shown in Fig. 2b. We can associate with it 
a stochastic-like variable 

:F{t) = -^{Q{t),P{t),x{t)), (17) 

which for our linearized Hamiltonian is simply the perturbation term 

T = — yy = —QiQl- It can be interpreted as the generalized force that acts on 
the boundary of the 2D well. It may have a non-zero average ("conservative" part) 
but below we are interested only in its fluctuations. 

In order to characterize the fluctuations of J^{t) we introduce the autocorrelation 
function C(t) 

C(r) = (^(t)^(t + r))-(^2^. (18) 

The angular brackets denote an averaging which is either micro-canonical over 
some initial conditions (Q(0), -P(O)) or temporal (due to the assumed ergodicity). 
The power spectrum for the 2D well model is shown in Fig.4 (see solid line). 
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Fig. 4. The classical power-spectrum of the model (14). The classical cut-off frequency 
Wci — 7 is indicated by perpendicular dashed lines. 

For generic chaotic systems (described by smooth Hamiltonians), the fluctuations 
are characterized by a short correlation time r^, after which the correlations are neg- 
ligible. In generic circumstances r^i is essentially the ergodic time. For our model 
system Td ~ 1. 

The power spectrum of the fluctuations C (to) is defined by a Fourier transform: 

C{ij) = / C(r) exp(2cur)dr . (19) 



This power spectrum is characterized by a cut-off frequency Ud which is inverse 
proportional to the classical correlation time 



27r 



(20) 



Indeed in the case of our model system we get Uci ~ 7 which is in agreement with 
Fig.4. 

The implication of having a short but non- vanishing classical correlation time Td is 
having large but finite bandwidth in the perturbation matrix B. This follows from 
the identity 

Ciu) = Y: \Bnrn\'2n6 L - EllZ^) , (21) 



which implies 



( I Bjim I ) 



h 



A ~ / En — Em 

-C \ UJ = 



(22) 



271^1 \ h 

Hence the matrix elements of the perturbation matrix B are extremely small outside 
of a band of width h = hud/ A. 

In the inset of Figure 5 we show a snapshot of the perturbation matrix It 
clearly shows a band-structure. At the same figure we also display the band-profiles 
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Fig. 5. The band-profile (27rfi/A) • |-B„mP versus uj = {En — Ejn)/his compared with 
the classical power spectrum (7(0;). Inset: a snapshot of the perturbation matrix B. 

for different values of h. A good agreement with the classical power spectrum C (lu) 
is evident. 

It is important to realize that upon quantization we end up with two distinct energy 
scales. One is obviously the mean level spacing (see previous subsection) 



(23) 



where the dimensionality is c? = 2 in case of our model system. The other energy 
scale is the bandwidth 

Ah = = bA . (24) 

Tcl 

This energy scale is also known in the corresponding literature as the "non-universal" 
energy scale [28], or (in case of diffusive motion) as the Thouless energy [29] ^ . 
One has to notice that deep in the semiclassical limit /i — > these two energy scales 
differ enormously from one another (provided d > 2). We shall see in the follow- 
ing sections that this scale separation has dramatic consequences in the theory of 
driven quantum systems. 



4.3 Distribution of couplings 



We investigate further the statistical properties of the matrix elements Bnm of the 
perturbation matrix, by studying their distribution. RMT assumes that upon appro- 
priate unfolding they must be distributed in a Gaussian maimer. The 'unfolding' 
aims to remove system specific properties and reveal the underlying universality. 
It is done by normalizing the matrix elements with the local standard deviation 
a — a/(|B„^P) related through Eq. (22) with the classical power spectrum C{u}). 



^ The dimensionless parameter b scales like b on h and in the frame of mesoscopic 

systems is recognized as the dimensionless Thouless conductance [29]. 
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Fig. 6. Distribution of matrix elements q around E = 3 rescaled with the averaged 
band-profile . The solid black line corresponds to a Gaussian distribution with unit vari- 
ance while the dashed-dotted line corresponds to a fit from Eq. (25) with a fitting parameter 
A'^ = 7.8. The quantization corresponds to h = 0.03. 



The existing literature is not conclusive about the distribution of the normalized 

matrix elements q = Bnm/o-. Specifically, Berry [30] and more recently Prosen 
[31,32], claimed that V{q) should be Gaussian. On the other hand, Austin and 
Wilkinson [33] have found that the Gaussian is approached only in the limit of 
high quantum numbers while for small numbers, i.e., low energies, a different dis- 
tribution applies, namely 



r(^) / ^2X(^-3)/2 

^*^^^^7^Mx^l'"^J ■ ^^^^ 

This is the distribution of the elements of an A^-dimensional vector, distributed 
randomly over the surface of an iV-dimensional sphere of radius V^. For — > oo 
this distribution approaches a Gaussian. 

The distribution V{q) for our model is reported in Figure 6. The solid line corre- 
sponds to a Gaussian of unit variance while the dashed-dotted line is obtained by 
fitting Eq. (25) to the numerical data using iV as a fitting parameter. We observe that 
the Gaussian resembles better our numerical data although deviations, especially 
for matrix elements close to zero, can be clearly seen. We attribute these deviations 
to the existence of the tiny stability islands in the phase space. Trajectories started 
in those islands cannot reach the chaotic sea and vice versa. Quantum mechanically 
the consequence of this would be vanishing matrix elements B^m which represent 
the classically forbidden transitions. 
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4. 4 RMT modeling 



It was the idea of Wigner [12,13] more than forty years ago, to study a simplified 
model, where the Hamiltonian is given by Eq. (8), and where B is a Banded Ran- 
dom Matrix (BRM) [34,35,36]. The diagonal matrix E has elements which are the 
ordered energies {E^}, with mean level spacing A. The perturbation matrix B has 
a rectangular band-profile of band-size h. Within the band < \n ~ m\ < b the 
elements are independent random variables given by a Gaussian distribution with 
zero mean and a variance cr^ = Outside the band they vanish. We refer 

to this model as the Wigner BRM model (WBRM). 

Given the band-profile, we can use Eq.(22) in reverse direction to calculate the 
correlation function C(r). For the WBRM model we get 

C(t) = 2(7^6 sine (r/rd) , (26) 

where — h/ A^. Thus, there are three parameters (A, b, a) that define the WBRM 
model. 

The WBRM model can be regarded as a simplified local description of a true Hamil- 
tonian matrix. This approach is attractive both analytically and numerically. Ana- 
lytical calculations are greatly simplified by the assumption that the off-diagonal 
terms can be treated as independent random numbers. Also from a numerical point 
of view it is quite a tough task to calculate the true matrix elements of the B matrix. 
It requires a preliminary step where the chaotic Ho is diagonalized. Due to memory 
limitations one ends up with quite small matrices. For the PuUen-Edmonds model 
we were able to handle matrices of final size N — 4000 maximum. This should 
be contrasted with the WBRM simulations, where using self -expanding algorithm 
[37,17] we were able to handle system sizes upto N — 100000 along with signifi- 
cantly reduced CPU time. 

We would like to stress again that the underlying assumption of WBRM, namely 
that the off-diagonal elements are uncorrelated random numbers, has to be treated 
with extreme care. 

The WBRM model involves an additional simplification. Namely, one assumes that 
B has a rectangular band-profile. A simple inspection of the band-profile of our 
model Eq. (14) shows that this is not the case (see Fig. 5). We eliminate this simpli- 
fication by introducing a RMT model that is even closer to the dynamical one. To 
this end, we randomize the signs of the off-diagonal elements of the perturbation 
matrix B keeping its band-structure intact. This procedure leads to a random model 
that exhibits only universal properties while it lack any semiclassical limit. We will 
refer to it as the effective banded random matrix model (EBRM). 



15 



5 The Parametric Evolution of the Eigenfunctions 



As we change the parameter 6x in the Hamiltonian Eq. (8), the instantaneous eigen- 
states {\n{x))} evolve and undergo structural changes. In order to understand the 
actual dynamics, it is important to understand these structural changes. This leads 
to the introduction of 

P{n\m) ^\{n{x)\m{xQ))\^ , (27) 

which is easier to analyze than P((n|no). Up to some trivial scaling and shifting 
P{n\m) is essentially the local density of states (LDoS): 

P{E\m) = Y,\{n{x)\m{xo))ME - E^) . (28) 

n 

The averaged distribution P(r) is defined in complete analogy with the definition of 
Pt{r). Namely, we use the notation r = n — m, and average over several m states 
with roughly the same energy E^ ~ E. 

Generically P(r) undergoes the following structural changes as a function of grow- 
ing Sx. We first summarize the generic picture, which involves the parametric scales 
Ec and Sprt. and the approximations Pfopt, Ppn, and P^c- Then we discuss how to de- 
termine these scales, and what these approximations are. 

• The first order perturbative theory (FOPT) regime is defined as the range 6x < Sc 
where we can use FOPT to get an approximation that we denote as P() ^ Pfopt- 

• The (extended) perturbative regime is defined as the range ec< Sx < e^^ where 
we can use perturbation theory (to an infinite order) to get a meaningful approx- 
imation that we denote as P() a; Ppn. Obviously Ppn reduces to Pfopt in the 
FOPT regime. 

• The non-perturbative regime is defined as the range 5x > £prt where perturbation 
theory becomes non-applicable. In this regime we have to use either RMT or 
semiclassics in order to get an approximation that we denote as P() fa Psc. 

Irrespective of these structural changes, it can be proved that the variance of P(r) 
is strictly linear and given by the expression 

5E{5x) = ^/Ofi) 6x = 6Ed . (29) 

The only assumption that underlines this statement is 6x -C ^ci- It reflects the linear 
departure of the energy surfaces. 
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5.1 Approximations for P{n\m) 



The simplest regime is obviously the FOPT regime where, for P(n|m), we can use 
the standard textbook approximation PYiopY{n\m) ^ 1 for n — m, while 

PFOPT(n|m) = — , (30) 

for n 7^ m. If outside of the band we have Bnm = 0, as in the WBRM model, then 
-Pfopt(^) = for |r| > 6. To find the higher order tails (outside of the band) we 
have to go to higher orders in perturbation theory. Obviously this approximation 
makes sense only as long as 6x < where 

£e = A/(7 ~ ^ (31) 

and d is the degrees of freedom of our system {d — 2 for the 2D well model). 

If bx > Ec but not too large then we still have tail regions which are described by 
FOPT This is a non-trivial observation which can be justified by using perturbation 
theory to infinite order. Then we can argue that a reasonable approximation is 

where F is evaluated by imposing normalization of Pprt (n | m) . In the case of WBRM 
model F = {a6x/AY x A. The appearance of F in the above expression cannot be 
obtained from any ^nzYe-orJer perturbation theory: Formally it requires summation 
to infinite order. Outside of the bandwidth the tails decay faster than exponentially. 
Note that Pprt(n|m) is a Lorentzian in the case of a flat bandwidth (WBRM model), 
while in the general case it can be described as a "core-tail" structure. 

Obviously the above approximation makes sense only as long as T{6x) < A5. This 
expression assumes that the bandwidth A5 is sharply defined, as in the WBRM 
model. By elimination this leads to the determination of epn, which in case of the 
WBRM model is simply 



h 



prt 



^ VbSc^ : . (33) 



Tel 



In more general cases the bandwidth is not sharply defined. Then we have to define 
the perturbative regime using a practical numerical procedure. The natural defini- 
tion that we adopt is as follows. We calculate the spreading SE{6x), which is a 
linear function. Then we calculate 5Ep^{5x), using Eq.(32)). This quantity always 
saturates for large 6x because of having finite bandwidth. We compare it to the 
exact 5E{5x), and define for instance as the 80% departure point. 
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Fig. 7. The parametric evolution of eigenstates of a WBRM model with a = I and b = 50: 
(a) Standard perturbative regime corresponding to e = 0.01, (b) Extended perturbative 
regime with e = 2 (c) Non-perturbative (ergodic) regime with e = 12 and (d) localized 
regime with e = 1. In (a-c) the mean level spacing A = 1 while in (d) A = 10^'^. The 
bandwidth Af, = A x 6 is indicated in all cases. In (b) the blue dashed line corresponds to 
a Lorentzian with F « 16 ^ A;, while in (a) we have F « 10~^ <^ A which therefore 
reduces to the standard FOPT result. 
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Fig. 8. The quantal profile P{r) for the 2DW model is compared with Ppn{r) and with 
the corresponding Prmt{i") of the EBRM model. The perturbation strength 6x is in (a) 
6x = 0.00035, (b) 6x = 0.04945 and (c) 6x = 0.29. We are using here the h = 0.012 
output. In the lower plot the classical LDoS profile Pci(^) is represented by a green heavy 
dashed line. 

What happens if perturbation theory completely fails? In the WBRM model the 
LDoS becomes semicircle: 



while in systems that have a semiclassical limit we expect to get 

Pso{n\m) = J^^^p4Q,P)p^(Q,P), (35) 

where pm{Q,P) and pn{Q,P) are the Wigner functions that correspond to the 
eigenstates \m{xo)) and \n{x)) respectively. 
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5.2 The P{n\m) in practice 



There are some findings tliat go beyond tlie above generic picture and, for com- 
pleteness, we mention them. The first one is the "localization regime" which is 
found in the case of the WBRM model for e > sioc- where 

£ioc = . (36) 

In this regime it is important to distinguish between the non-averaged P{n\m) and 
the averaged P(r) because the eigenfunctions are non-ergodic but rather localized. 
This localization is not reflected in the LDoS which is still a semicircle. A typi- 
cal eigenstate is exponentially localized within an energy range 5E^ — much 
smaller than 5E^\. The localization length is { ~ 6^. In actual physical applications 
it is not clear whether there is such a type of localization. The above scenario for 
the WBRM model is summarized in Fig. 7 where we plot P{n\m) in the various 
regimes. The localized regime is not an issue in the present work and therefore we 
will no further be concerned with it. 

The other deviation from the generic scenario, is the appearance of a non-universal 

"twilight regime" which can be found for some quantized systems [38]. In this 
regime a co-existence of a perturbative and a semiclassical structure can be ob- 
served. For the PuUen-Edmonds model (14) there is no such distinct regime. 

For the Hamiltonian model described by Eq. (14) the borders between the regimes 
can be estimated [19]. Namely Sc ~ and £prt ~ 5.3/1. In Fig. 8 we report 

the parametric evolution of the eigenstates for the Hamiltonian model of Eqs. (14) 
and we compare the outcomes with the results of the EBRM model [19]. Despite 
the overall quantitative agreement, some differences can be detected: 

• In the FOPT regime (see Fig. 8a), the RMT strategy fails in the far tails regime 
A X |r| > Afe where system specific interference phenomena become important. 

• In the extended perturbative regime (see Fig. 8b) the line-shape of the averaged 
wavefunction P{n\m) is different from Lorentzian. Still the general features of 
Pprt (core-tail structure) can be detected. In a sense, Wigner's Lorentzian (32) is a 
special case of core-tail structure. Finally, as in the standard perturbative regime 
one observes that the far-tails are dominated by either destructive interference 
(left tail), or by constructive interference (right tail). 

• Deep in the non-perturbative regime i s > s^n ) the overlaps P(n|m) are well 
approximated by the semiclassical expression. The exact shape is determined by 
simple classical considerations [19,39]. This is in contrast to the WBRM model 
which does not have a classical limit. 
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6 Linear Response Theory 



The definition of regimes for driven systems is more complicated than the corre- 
sponding definition in case of LDoS theory. It is clear that for short times we al- 
ways can use time-dependent FOPT. The question is, of course, what happens next. 
There we have to distinguish between two types of scenarios. One type of scenario 
is wavepacket dynamics for which the dynamics is a transient from a preparation 
state to some new ergodic state. The second type of scenario is persistent driving, 
either linear driving {x — e) or periodic driving {x{t) — e sin(Qt)). In the latter 
case the strength of the perturbation depends on the rate of the driving, not just 
on the amplitude. The relevant question is whether the long time dynamics can 
be deduced from the short time analysis. To say that the dynamics is of perturba- 
tive nature means that the short time dynamics can be deduced from FOPT, while 
the long time dynamics can be deduced on the basis of a Markovian (stochastic) 
assumption. The best known example is the derivation of the exponential Wigner 
law for the decay of metastable state. The Fermi-Golden-Rule (FOR) is used to 
determine the initial rate for the escaping process, and then the long-time result is 
extrapolated by assuming that the decay proceeds in a stochastic-like manner. Simi- 
lar reasoning is used in deriving the Pauli master equation which is used to describe 
the stochastic-like transitions between the energy levels in atomic systems. 

A related question to the issue of regimes is the validity of Linear Response Theory 
(LRT). In order to avoid ambiguities we adopt here a practical definition. When- 
ever the result of the calculation depends only on the two point correlation func- 
tion C(r), or equivalently only on the band-profile of the perturbation (which is 
described by C{uj)), then we refer to it as "LRT". This implies that higher order 
correlations are not expressed. There is a (wrong) tendency to associate LRT with 
FOPT. In fact the validity of LRT is not simply related to FOPT. We shall clarify 
this issue in the next section. 

For both 5E{t) and V{t) we have "LRT formulas" which we discuss in the next 
sections. Writing the driving pulse as 5x{t) = ef{t) for the spreading we get: 

/CO rji.] ^ _ 
—C{u)F,{u) , (37) 

while for the survival probability we have 

Two spectral functions are involved: One is the power spectrum C {u) of the fluctu- 
ations defined in Eq. (19), and the other Ft{uj) is the spectral content of the driving 
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pulse which is defined as 



[' dt'f(t')e-'^'' 
Jo 



(39) 



Here we summarize the main observations regarding the nature of wavepacket dy- 
namics in the various regimes: 

• FOPT regime: In this regime V{t) ~ 1 for all time, indicating that all probability 

is all the time concentrated on the initial level. An alternative way to identify this 
regime is from (5£'core(^) which is trivially equal to A. 

• Extended perturbative regime: The appearance of a core-tail structure which is 
characterized by separation of scales A <^ 5Ecoie.{t) -C 5E{t) <C A^. The core 
is of non-perturbative nature, but the variance SE'^(t) is still dominated by the 
tails. The latter are described by perturbation theory. 

• Non-perturbative regime: The existence of this regime is associated with having 
the finite energy scale A^. It is characterized by A;, <^ 6Ecoreit) ~ 5E{t). As 
implied by the terminology, perturbation theory (to any order) is not a valid tool 
for the analysis of the energy spreading. Note that in this regime, the spreading 
profile is characterized by a single energy scale {5E ~ SEcor^). 



6.1 The energy spreading SE{t) 



Of special importance for understanding quantum dissipation is the theory for the 
variance 5E'^{t) of the energy spreading. Having SE{t) oc e means linear response. 
If SE{t)/e depends on e, we call it "non-Unear response". In this paragraph we 
explain that linear response theory (LRT) is based on the "LRT formula" Eq.(37) 
for the spreading. This formula has a simple classical derivation (see Subsection 
6.2 below). 

From now on it goes without saying that we assume the classical conditions for the 
validity of Eq.(37) are satisfied (no h involved in such conditions). The question 
is what happens to the validity of LRT once we "quantize" the system. In previous 
publications[8, 10,1 1,19], we were able to argue the following: 

(A) The LRT formula can be trusted in the perturbative regime, with the exclusion 
of the adiabatic regime. 

(B) In the sudden limit the LRT formula can also be trusted in the non-perturbative 

regime. 

(C) In general the LRT formula cannot be trusted in the non-perturbative regime. 

(D) The LRT formula can be trusted deep in the non-perturbative regime, provided 
the system has a classical limit. 

For a system that does not have a classical limit (Wigner model) we were able to 
demonstrate [8,10,11] that LRT fails in the non-perturbative regime. Namely, for 



22 



the WBRM model the response SE{t) / e becomes e dependent for large e, meaning 
that the response is non-linear. Hence the statement in item (C) above has been 
established. We had argued that the observed non-linear response is the result of a 
quantal non-perturbative effect. Do we have a similar type of non-linear response 
in the case of quantized chaotic systems? The statement in item (D) above seems to 
suggest that the observation of such non-linearity is not likely. Still, it was argued in 
[11] that this does not exclude the possibility of observing a "weak" non-linearity. 

The immediate (naive) tendency is to regard LRT as the outcome of quantum me- 
chanical first order perturbation theory (FOPT). In fact the regimes of validity of 
FOPT and of LRT do not coincide. On the one hand we have the adiabatic regime 
where FOPT is valid as a leading order description, but not for response calculation. 
On the other hand, the validity of Eq.(37) goes well beyond FOPT. This leads to 
the (correct) identification [7,8,11] of what we call the "perturbative regime". The 
border of this regime is determined by the energy scale A;,, while A is not involved. 
Outside of the perturbative regime we cannot trust the LRT formula. However, as 
we further explain below, the fact that Eq.(37) is not valid in the non-perturbative 
regime, does not imply that it fails there. 

We stress again that one should distinguish between "non-perturbative response" 
and "non-linear response". These are not synonyms. As we explain in the next 
paragraph, the adiabatic regime is "perturbative" but "non-linear", while the semi- 
classical limit is "non-perturbative" but "linear". 

In the adiabatic regime, FOPT implies zero probability to make a transitions to 
other levels. Therefore, to the extent that we can trust the adiabatic approxima- 
tion, all probability remains concentrated on the initial level. Thus, in the adiabatic 
regime, Eq.(37) is not a valid formula: It is essential to use higher orders of per- 
turbation theory, and possibly non-perturbative corrections (Landau-Zener [1,2]), 
in order to calculate the response. Still, FOPT provides a meaningful leading order 
description of the dynamics (i.e. having no transitions), and therefore we do not 
regard the adiabatic non-linear regime as "non-perturbative". 

In the non-perturbative regime the evolution of Pt{n\m) cannot be extracted from 
perturbation theory: not in leading order, neither in any order. Still it does not neces- 
sarily imply a non-linear response. On the contrary: The semiclassical limit is con- 
tained in the deep non-perturbative regime [8,1 1]. There, the LRT formula Eq.(37) 
is in fact valid. But its validity is not a consequence of perturbation theory, but 
rather the consequence of quantal-classical correspondence (QCC). 

In the next subsection we will present a classical derivation of the general LRT 
expression (37). In Subsection 6.3 we derive it using first order perturbation theory 
(FOPT). In Subsection 6.5 we derive the corresponding FOPT expression for the 
survival probability. 
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6.2 Classical LRT derivation for SE{t) 



The classical evolution of E{t) — T-C{Q{t), P{t)) can be derived from Hamilton 
equations. Namely, 

^ = [n, nu + ^ = -ef{t)m , (40) 

where [-jpB indicates the Poisson Brackets. Integration of Eq. (40) leads to 

E(t) - E(0) = -6 /* J^(t')jXt')dt' . (41) 
Jo 

Taking a micro-canonical average over initial conditions we obtain the following 
expression for the variance 

SE\t) = t Cit' - t") f(t')f(t")dt'dt" , (42) 
Jo 

which can be re-written in the form of (37). 

One extreme special case of Eq.(37) is the sudden limit for which /(t) is a step 
function. Such evolution is equivalent to the LDoS studies of Section 5. In this case 
Ft{uj) = 1, and accordingly 



5Ea = £ X ^JC{0) ["sudden" case] . (43) 



Another extreme special case is the response for persistent (either linear or periodic) 
driving of a system with an extremely short correlation time. In such case Ft{uj) 
becomes a narrow function with a weight that grows linearly in time. For linear 
driving (f{t) = t) we get Ft{uj) = t x 27r5{uj). This implies diffusive behavior: 



5E{t) = ^J2DEt ["Kubo" case] , (44) 

where De oc is the diffusion coefficient. The expression for De as an inte- 
gral over the correlation function is known in the corresponding literature either as 
Kubo formula, or as Einstein relation, and is the comer stone of the Fluctuation- 
Dissipation relation. 



6.3 Quantum LRT derivation for 5E{t) 



The quantum mechanical derivation looks like an exercise in first order perturbation 
theory. In fact a proper derivation that extends and clarifies the regime where the 
result is applicable requires infinite order. If we want to keep a complete analogy 



24 



with the classical derivation we should work in the adiabatic basis [7]. (For a brief 
derivation see Appendix D of [9]). 



In the following presentation we work in a "fixed basis" and assume f{t) — f{0) — 
0. We use the standard textbook FOPT expression for the transition probability 
from an initial state m to any other state n. This is followed by integration by parts. 
Namely, 



Pt{n\m) = —\B, 



2 







_£:iR ,2: 



(45) 



where ujnm = {En — Em) I h. Now we calculate the variance and use Eq. (22) so as 
to get 



SE\t) = J2Pt{n\m){En-Emf 



(46) 



6.4 Restricted QCC 



The FOPT result for SE{t) is exactly the same as the classical expression Eq. (37). 

It is important to realize that there is no /i-dependence in the above formula. This 
correspondence does not hold for the higher A:— moments of the energy distribution. 
If we use the above FOPT procedure we get that the latter scale as h'^~'^. 

We call the quantum-classical correspondence for the second moment "restricted 
QCC". It is a very robust correspondence [11]. This should be contrasted with "de- 
tailed QCC" that applies only in the semiclassical regime where Pt{n\m) can be 
approximated by a classical result (and not by a perturbative result). 



6. 5 Quantum LET derivation for V (t) 



With the validity of FOPT assumed we can also calculate the time-decay of the 
survival probability V{t). From Eq. (45) we get: 

E P,(n|»)^.7_;|cMg). (47) 
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Assuming that V{t) = 1 —pit) can be extrapolated in a "stochastic" fashion we get 
Eq. (38). Another way to write the final formula is as follows: 



V{t) = exp 



^ f f Cit'-t")6xit')6xit")dt'dt 

Jo Jo 



(48) 



For constant perturbation (wavepacket dynamics) and assuming long times we ob- 
tain the Wigner decay. 



V{t) = exp 



^\ Ciu^Q) X t 
hj 



(49) 



which can be regarded as a special case of Fermi-Golden-Rule. 

6.6 Note on V{t) for a time reversal scenario 

The "LRT formula" for V{t) in the case of "driving reversal scenario" is 

2 pT fT 



Pdr(^) = exp 



hJ Jo Jo 



C{t'-t")f{t')f{t")dt'dt" 



(50) 



where we assumed the simplest scenario with f{t) = 1 for < t < (T/2) and 
f{t) = —1 for (T/2) < t < T. It is interesting to make a comparison with the anal- 
ogous result in case of "time reversal scenario". 

The well known Feynman- Vernon influence functional has the following approxi- 
mation: 



F[xA,XB]^{^\U[xB]-^U[xAm 
1 f* /■* 



(51) 



= exp 



2n? Jo Jo 



C{t'-t"){xB{t')-XA{t")fdt'dt" 



This expression is in fact exact in the case of harmonic bath, and assuming thermal 
averaging over the initial state. Otherwise it should be regarded as an extrapolated 
version of leading order perturbation theory (as obtained in the interaction picture). 
What people call nowadays "fidelity" or "Loschmidt echo" is in fact a special case 
of the above expression which is defined by setting t — T/2 and xa — s/2 while 
xb — —e/2. Thus 



VTK{t) = \F[xA,XB\? 

-P\2 fT/2 fT/2 

M / c{t'-t")dt'dt" 

Jo Jo 



(52) 



: exp 
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Assuming a very short correlation time one obtains 




) 



2 



Vtr{T) = exp 



(53) 



which again can be regarded as a special variation of the Fermi-Golden-Rule (but 
note the pre-factor 1/2). 

6. 7 The survival probability and the LDoS 

For constant perturbation it is useful to remember that V{t) LDoS as follows: 



This implies that a Wigner decay is associated with a Lorentzian approximation 
for the LDoS. In the non-perturbative regime the LDoS is not a Lorentzian, and 
therefore one should not expect an exponential. In the semiclassical regime the 
LDoS shows system specific features and therefore the decay of V{t) becomes 
non-universal. 



7 Wavepacket Dynamics for Constant Perturbation 

The first evolution scheme that we are investigating here is the so-called wavepacket 
dynamics. The classical picture is quite clear [17,18]: The initial preparation is 
assumed to be a micro-canonical distribution that is supported by the energy surface 
T-Co{Q,P) — E{0). Taking to be a generator for the classical dynamics, the 
phase-space distribution spreads away from the initial surface for t > 0. 'Points' 
of the evolving distribution move upon the energy surfaces of ^.{Q, P). Thus, the 
energy E{t) = Ho{Q{t), P{t)) of the evolving distributions spreads with time. 
Using the LRT formula Eq.(39) for rectangular pulse f{t') — 1 for < t' < t we 
get 



7^(i) = |(n(a;o)|e-*^(^)*/VW)f 

= ^e-'-^-(^)*/'^|(m(a;)|n(a;o))|' 



2 



m 




(54) 




(55) 



and hence 




(56) 
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Fig. 9. Simulations of wavepacket dynamics for the 2DW model (left panel) and for the 
corresponding EBRM model (right panel). The energy spreading SE{t) (normalized with 

respect to the perturbation strength e) is plotted as a function of time for various perturba- 
tion strengths e corresponding to different hne-types (the same in both panels). The classical 
spreading 5Ec\ (t) (thick dashed line) is plotted in both panels as a reference. 

For short times t <^ we can expand the correlation function as C{t) ^ C(0) — 
|C"(0)t^, leading to a ballistic evolution. Then, for t 3> r^, due to ergodicity, a 
'steady-state distribution' appears, where the evolving 'points' occupy an 'energy 
shell' in phase-space. The thickness of this energy shell equals SEa. Thus, we have 
a crossover from ballistic energy spreading to saturation: 



Figure 9 shows the classical energy spreading (heavy dashed line) for the 2DW 
model. In agreement with Eq. (57) we see that SEc:\{t) is first ballistic and then sat- 
urates. The classical dynamics is fully characterized by the two classical parameters 
Tci and SEa. 

7.1 The quantum dynamics 

Let us now look at the quantized 2DW model. The quantum mechanical data are 
reported in Fig. 9 (left panel) where different curves correspond to various pertur- 
bation strengths e. As in the classical case (heavy dashed-line) we observe an initial 
ballistic-like spreading [18] followed by saturation. This could lead to the wrong 
impression that the classical and the quantum spreading are of the same nature. 
However, this is definitely not the case. 

In order to detect the different nature of quantum ballistic-like spreading, one has 
to inquire measures that are sensitive to the structure of the profile, such as the 
core- width 6Ecore{t). In Fig. 10 we present our numerical data for the 2DW model. 
If the spreading were of a classical type, it would imply that the spreading pro- 




(57) 
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Fig. 10. Simulations of wavepacket dynamics for the 2DW model. The evolution of the 
(normalized) core width SEcore {t) is plotted as a function of time. The classical expectation 
is represented by a thick dashed Une for the sake of comparison. As e becomes larger it is 
approached more and more. We use the same set of parameters as in Fig. 9. 

file is characterized by a single energy scale. In such a case we would expect that 
SEcore{t) ~ 6E{t). Indeed this is the case for e > £prt with the exclusion of very 
short times: The larger e is the shorter the quantal transient becomes. In the per- 
turbative regimes, in contrast to the semiclassical regime, we have a separation of 
energy scales 6Ecoreit) <^ 5E{t). In the perturbative regimes 5E{t) is determined 
by the tails, and it is not sensitive to the size of the 'core' region. 

Using the LRT formula for V{t) we get, for short times {t <^ Td) during the 
ballistic-like stage 



while for long times {t » Td) we have the FGR decay of Eq.(49). Can we trust 
these expressions? Obviously FOPT can be trusted as long as V{t) ~ 1. This can 
be converted into an inequality t < tp^ where 



The power u = 1 applies to the non-perturbative regime where the breakdown of 
V{t) happens to be before r^. The power u = 2 applies to the perturbative regime 
where the breakdown of V{t) happens after at tp^t = h/T, i.e. after the ballistic- 
like stage. 

The long term behavior of V{t) in the non-perturbative regime is not the Wigner 
decay. It can be obtained by Fourier transform of the LDoS. In the non-perturbative 
regime the LDoS is characterized by the single energy scale SE^i oc 5x. Hence the 
decay in this regime is characterized by a semiclassical time scale 27rh/ 5Ec\. 




(58) 




(59) 
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7.2 The EBRM dynamics 



Next we investigate the applicability of the RMT approach to describe wavepacket 
dynamics [17,18] and specifically the energy spreading SE{t). At first glance, we 
might be tempted to speculate that RMT should be able, at least as far as ()E{t) is 
concerned, to describe the actual quantum picture. After all, we have seen in Sub- 
section 6. 1 that the quantum mechanical LRT formula (46) for the energy spreading 
involves as its only input the classical power spectrum C{u;). Thus we would ex- 
pect that an effective RMT model with the same band-profile would lead to the 
same SE{t). 

However, things are not so trivial. In Figure 9 we show the numerical results for 
the EBRM model ^ . In the standard and in the extended perturbative regimes we 
observe a good agreement with Eq.(46). This is not surprising as the theoretical 
prediction was derived via FOPT, where correlations between off-diagonal ele- 
ments are not important. In this sense the equivalence of the 2DW model and 
the EBRM model is trivial in these regimes. But as soon as we enter the non- 
perturbative regime, the spreading 6E{t) shows a qualitatively different behavior 
from the one predicted by LRT: After an initial ballistic spreading, we observe a 
premature crossover to a diffusive behavior 

SE{t) = ^2D^t . (60) 

The origin of the diffusive behavior can be understood in the following way. Up 
to time iprt the spreading 5E{t) is described accurately by the FOPT result (46). 

At i ~ iprt the evolving distribution becomes as wide as the bandwidth, and we 
have 5Ecove ^ 6E Ab rather than SEcore 6E <^ At,. We recall that in the 
non-perturbative regime FOPT is subjected to a breakdown before reaching satura- 
tion. The following simple heuristic picture turns out to be correct. Namely, once 
the mechanism for ballistic-like spreading disappears, a stochastic-like behavior 
takes its place. The stochastic energy spreading is similar to a random-walk pro- 
cess where the step size is of the order A;,, with transient time tpn. Therefore we 
have a diffusive behavior 5E{ty = 2D^t with 

De = C-Al/tpn ■ AH^/^ea/h oc H (61) 

where C is some numerical pre-factor. This diffusion is not of classical nature, 
since in the — > limit we get De — > 0. The diffusion can go on until the en- 
ergy spreading profile ergodically covers the whole energy shell and saturates to a 
classical-like steady state distribution. The time terg for which we get ergodization 
is characterized by the condition (De^)^^^ < ^-^ci, leading to 

i,,g = h-^l'^hea/A^ (X 1/h. (62) 
^ The same qualitative results were found also for the prototype WBRM model, see [17]. 
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Fig. 11. A diagram that illustrates the various time scales in wavepacket dynamics, de- 
pending on the strength of the perturbation e. The diagram on the left refers to the WBRM 
model, while that on the right is for a quantized system that has a classical limit. The two 
cases differ in the non-perturbative regime (large e): In the case of a quantized model we 
have a genuine ballistic behavior which reflects detailed QCC, while in the RMT case we 
have a diffusive stage. In the latter case the times scale Udn marks the crossover from re- 
versible to non-reversible diffusion. This time scale can be detected in a driving reversal 
scenario as explained in the next section. For further discussion of this diagram see the 
text, and in particular the concluding section of this paper. 

For completeness we note that for e > eioc there is no ergodization but rather 
dynamical ("Anderson" type) localization. Hence, in the latter case, terg is replaced 
by the break-time t\,±. The various regimes and time scales are illustrated by the 
diagram presented in Fig. 1 1 . 



8 Driving Reversal Scenario 



A thorough understanding of the one-period driving reversal scenario [10] is both 
important within itself, and for constituting a bridge towards a theory dealing with 
the response to periodic driving [8]. In the following subsection we present our 
results for the prototype WBRM model, while in Subsection 8.2 we consider the 
2DW model and compare it to the corresponding EBRM model. The EBRM is 
better for the purpose of making comparisons with the 2DW, while the WBRM 
is better for the sake of quantitative analysis (the "physics" of the EBRM and the 
WBRM models is, of course, the same). 

The quantities that monopolize our interest are the energy spreading SE{t) and the 
survival probability P (t) . In Figs . 1 2 and 1 3 we present representative plots . From a 
large collection of such data that collectively span a very wide range of parameters, 
we extract results for 6E(T), for V{T), and for the corresponding compensation 
times. These are presented in Figs. 12,13,14,15,16 and 17. 
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Fig. 12. Simulations of driving reversal for the 2DW model (left panels) and for the corre- 
sponding EBRM model (right panels). In the upper row the (normalized) energy spreading 
6E{t) is plotted as a function of time for representative values e while T = 0.48. In the 
lower row the (normalized) energy spreading 5E{T) at the end of the cycle is plotted versus 
T for representative values of e. 

8. 1 Driving Reversal Scenario: RMT Case 



8.1.1 LRT for the energy spreading 

Assuming that the driving reversal happens att = T/2, the spectral content Ft{uj) 
for T/2 <t<T is 



1 - 2e~*'^^/^ + e**^* 



(63) 



Inserting Eq. (63) into Eq. (37) we get 



5E{t) = e X J6C(0) + 2C{t) - 4C(|) - 4C(t-|: 



(64) 



For the WBRM model we can substitute in Eq. (64) the exact expression Eq. (26) 
for the correlation function, and get 



5E{t) = 2ea x 



\ 



t T t—-. 

3b + 6sinc( — ) — 2teinc( ) — 2teinc(- 

^ci 2r ci 



(65) 
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Fig. 13. Simulations of driving reversal for the 2DW model (left panels) and for the cor- 
responding EBRM model (right panels). In the upper row the (normalized) transition prob- 
ability p{t) = 1 — V{t) is plotted as a function of time for representative values e. The 
period of the driving is T = 0.48. In the lower row the survival probability V{T) at the end 
of the cycle is plotted versus T for representative values of e. In all cases we are using the 
same symbols as in the upper left panel of Fig. 12. 

We can also find the compensation time tf by minimizing Eq. (64) with respect to 
t. For the WBRM model we have 



2 cos 


"T/2-t" 




L Tel -1 



+ 



2 sin 


'T/2-t' 







cos 



+ 



t 



rei(T/2 - t) (T/2 - t) 
which can be solved numerically to get tf. 
The spreading width at the end of the period is 



tTcl 



t 

Tcl 



(66) 



SE{T) = ex W6C(0) + 2C{T) - 8C(^)) . (67) 



It is important to realize that the dimensional parameters in this LRT analysis are 
determined by the time scale and by the energy scale SEd. This means that we 
have a scaling relation (using units such that a = A = h = 1) 



Kkt m . (68) 
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Deviation from this scaling relation implies a non-perturbative effect that goes be- 
yond LRT. 

The LRT scaling is verified nicely by our numerical data (see upper panels of 
Fig. 14 ). The values of perturbation strength for which the LRT results are appli- 
cable correspond to s < Sp^. In the same Figure we also plot the whole analytical 
expression (67) for the spreading 5E(T). Similarly in Fig. 14 (lower panels) we 
present our results for the compensation time tf. All the data fall one on top of the 
other once we rescale them. It is important to realize that the LRT scaling relation 
implies that the compensation time is independent of the perturbation strength e. 
It is determined only by the classical correlation time Tc\. In the same figure we also 
present the resulting analytical result (heavy-dashed line) which had been obtained 
via Eq.(66). An excellent agreement with our data is evident. 



8.1.2 Energy spreading in the non-perturbative regime 

We turn now to discuss the dynamics in the non-perturbative regime, which is our 
main interest. In the absence of driving reversal (see Subsection 7.2) we obtain 
diffusion (SE{t) oc y^) for t > tpn, where 



If (T/2) < tprt, this non-perturbative diffusion does not have a chance to develop, 
and therefore we can still trust Eq. (64). So the interesting case is {T/2) > tpn, 
which means large enough e. In the following analysis we distinguish between two 
stages in the non-perturbative diffusion process. The first stage (tp^ < t < tsdn) is 
reversible, while the second stage (t > tsdn) is irreversible. For much longer time 
scales we have recurrences or localization, which are not the issue of this paper. The 
new time scale (^sdn) did not appear in our "wavepacket dynamics" study, because 
it can be detected only by time driving reversal experiment. 

The determination of the time scale t^dn is as follows. The diffusion coefficient is 
— A%^''^ae/h up to a numerical pre-factor. The diffusion law is 5E'^{t) = D^t. 
The diffusion process is reversible as long as E does not affect the relative phases 
of the participating energy levels. This means that the condition for reversibility is 
{5E{t) X t)/h<^ 1. The latter inequality can be written ast < t^dn, where 



It is extremely important to realize that without reversing the driving, the presence 
or the absence of E in the Hamiltonian cannot be detected. It is only by driving 
reversal that we can easily determine (as in the upper panels of Fig. 12) whether the 
diffusion process is reversible or irreversible. 



^prt = ^/ (Vbae) 



(69) 




(70) 
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Fig. 14. Simulations of driving reversal for the WBRM model. In the upper row the (scaled) 
energy spreading SE{T) at the end of the cycle is plotted against the (scaled) period T. In 
the lower panels the compensation time is plotted against the (scaled) period T. The panels 
on the left are for e values within the perturbative regime, while the panels on the left are 
for the non-perturbative regime. For the sake of comparison we plot the LRT expectation 
for 6 = 1 as a heavy dashed hne. 



The dimensional parameters in this analysis are naturally the time scale tsdn and the 
resolved energy scale h/T. Therefore we expect to have instead of the LRT scaling, 
a different "non-perturbative" scaling relation. Namely, 5E{T)/{h/T) should be 
related by a scaling function to T/t^^. Equivalently (using units such that 
a — A — h — l)it can be written as 



SE{T) 

65/6£l/3 



"■nprt 



(71) 



Obviously the non-perturbative scaling with respect to e^/^ goes beyond any impli- 
cations of perturbation theory. It is well verified by our numerical data (see upper 
right panel of Fig. 14). The values of perturbation strength for which this scaling 
applies correspond to s > Eprt- The existence of the tsdn scaling can also be verified 
in the lower right panel of Fig. 14, where we show that /T is by a scaling function 
related to h^/^e^'^T. 
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8.1.3 Decay ofV{t) in the FOPT regime 



We can substitute Eq. (63) for the spectral content Ft{uj) of the driving into the 
LRT formula Eq. (47), and come out with the following expression for the survival 
probability at the end of the period t — T 

V{T) fti ex^{-e^T%^) . (72) 

This is a super-Gaussian decay, which is quite different from the standard Gaus- 
sian decay Eq. (58) or any other results on reversibility that appear in literature 
[20,21,22,23,24]. We have verified that this expression is valid in the FOPT regime. 
See Figure 15 a. 

For the WBRM model, we get the following expression for p{t) after substituting 
the spectral content Ft{uj) given from Eq. (63) 



pit) = ^ X r d J ~ ^^"""^"^^ + """^^^^ ~ + ^^^^^^^^ . (73) 



The corresponding compensation time can be found after minimizing the above 
expression (corresponding to the maximization ofV{t) = 1 — p{t)) with respect to 
time t. This results in the following equation 

si{u;at)^2si(^uja(t-^^^ , (74) 

which has to be solved numerically in order to evaluate t^. Above si(a;) = Jq 
Our numerical data are reported in Fig. 16 together with the theoretical prediction 
(74). 



8.1.4 Decay ofV{t) in the Wigner regime 

We now turn to discuss V{t) in the "Wigner regime". By this we mean Ec < £ < 
£prt. This distinction does not appear in the 5E{t) analysis. The time evolution of 
SE{t) is dominated by the tails of the distribution and does not affect the "core" 
region. Therefore 6E{t) also agreed with LRT outside of the FOPT regime in the 
whole (extended) perturbative regime. But this is not the case with Vit), which is 
mainly influenced by the "core" dynamics. As a result in the "Wigner regime" we 
get different behavior compared with the FOPT regime. 

We look at the survival probability V{T) at the end of the driving period. In the 
Wigner regime, instead of the LRT-implied super-Gaussian decay, we find a Wigner- 
like decay: 

P(T) «e-^(^)^, (75) 
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where F ?a e^/A. In Figure 15b we present our numerical results for various per- 
turbation strengths in this regime. A nice overlap is observed once we rescale the 
time axis as £^ x T. We would like to emphasize once more that both in the standard 
and in the extended perturbative regimes the scaling law involves the perturbation 
strength e. This should be contrasted with the LRT scaling of SE{t). 

What about the compensation time ? A reasonable assumption is that it will ex- 
hibit a different scaling in the FOPT regime and in the Wigner regime (as is the 
case of V{t)). Namely, in the FOPT regime we would expect "LRT scaling" with 
Tci, while in the Wigner regime we would expect "Wigner scaling" with tpn — 
The latter is of non-perturbative nature and reflects the "core" dynamics. To our 
surprise we find that this is not the case. Our numerical data presented in Fig. 16 
show beyond any doubt that the "LRT scaling" applies within the whole (extended) 
perturbative regime, as in the case of tf, thus not invoking the perturbation strength 
e. We see that the FOPT expression (74) for shown as a heavy-dashed line de- 
scribes the numerical findings. 

We conclude that the compensation time U is mainly related to the dynamics of the 
tails, and hence can be deduced from the LRT analysis. 

8.1.5 Decay ofVit) in the non-perturbative regime 

Let us now turn to the non-perturbative regime (see Fig. 15c). As in the case of the 
spreading kernel 6E{T), the decay of V{T) is no longer captured by perturbation 
theory. Instead, we observe the same non-universal scaling with respect to e^^^ x T 
as in the case of 6E{T). 



The reason is that in the non-perturbative regime the two energy scales F and A;,, 
which were responsible for the difference between V{T) and 5E(T), lose their 
meaning. As a consequence, the spreading process involves only one time scale 
and the behavior of both V{T) and 6E{T) becomes similar, leading to the same 
scaling behavior. 



8.2 Driving Reversal Scenario: 2DW Case 

In the representative simulations of the 2DW model in Fig. 12 (upper left panel) 
we see that the spreading 5E{t) for T = 0.48 and various perturbation strengths 
e follows the LRT predictions very well. Fig. 12 (lower left panel) shows that the 
agreement with the LRT is observed for any value of the period T. This stands in 
clear contrast to the EBRM model shown in Fig. 12 (right panels). 




(76) 
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Fig. 15. Simulations of driving reversal for the WBRM model. The survival probabiUty 
V{T) at the end of the pulse is plotted against T in the (a) FOPT regime, (b) Wigner regime 
and (c) non-perturbative regime. In (a) the thick dashed line indicates the super-Gaussian 
decay (72) wliile (b-c) indicates a Wigner exponential decay (75). Various symbols corre- 
spond to different (e,b) values such that the e < ec in (a); ec < e < eprt in (b) and e > Cprt 
in (c). 
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Fig. 16. The compensation time is extracted for the same simulations of the previous fig- 
ures. Note that for both the FOPT regime and the Wigner regime we have "LRT scaling". 
The thick dashed Une corresponds to the LRT prediction Eq. (74) for 6 = 1. Various sym- 
bols correspond to different (e,b) values such that the e < Cc in (a); Cc < e < eprt in (b) and 
e > eprt in (c). 
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The agreement with LRT in the non-perturbative regime, as in the case of wavepacket 
dynamics, reflects detailed QCC. We recall that "to get into the non-perturbative 
regime" and "to make h small" means the same. All our simulations are done in 
a regime where LRT can be trusted at the classical (=non-perturbative) limit. It is 
only for RMT models that we observe a breakdown of LRT in the non-perturbative 
regime. 

What about V{T)1 This quantity has no classical analogue. Therefore QCC con- 
siderations are not applicable. Also LDoS considerations cannot help here. The 
one-to-one correspondence between the LDoS and the survival probability applies 
to the simple wavepacket dynamics scenario (constant perturbation). 

It is practically impossible to make a quantitative analysis of V{T) in the case 
of a real model because the band-profile is very structured and there are severe 
numerical limitations. Rather, what we can easily do is to compare the 2DW with 
the corresponding EBRM. Any difference between the two constitutes an indication 
for a non-perturbative effect. Representative simulations are presented in Figure 13. 

In Figure 17 we show the dependence of the compensation time on T for the 
EBRM model. We see very nice scaling behavior that indicates that our numerics 
(as far as V{T) is concerned!) is limited to the perturbative regime. We emphasize 
again that the physics of V{t) is very different from the physics of 5E{t). Therefore, 
this finding by itself should not be regarded as very surprising. A sharp crossover to 
a non-perturbative behavior can be expected for a "sharp" band-profile only (which 
is the WBRM and not the EBRM- see Fig. 16). 

Now we switch from the EBRM model to the 2DW model. Do we see any deviation 
from LRT scaling? The answer from Fig. 17 is clearly yes, as reflected by the e 

dependence of the curve. The effect is small, but "it is there". It indicates that the 
"body" of the probability distribution, in the case of the 2DW dynamics, does not 
evolve the same way as in the EBRM case. Indeed we know that the main part of 
the distribution evolves faster (in a ballistic fashion rather than diffusively), and 
therefore we observe lower values of tf. 

Assuming that the decay of P(T) is given by the exponential law, we extract the 
corresponding decay rates 7. It should be clear that the fitting is done merely in 
order to extract a numeric measure for the behavior of the decay. We would not 
like to suggest that the decay looks strictly exponential. The results are reported 
in Figure 18. We find that for e < e^^t, the decay rate 7(e) oc e^, as expected by 
Wigner's theory, while for e > Sp^t we find that 7 oc This linear dependence on 
e is essentially the same as in the corresponding wavepacket dynamics scenario. 
There it is clearly associated with the width 6Eci oc £ of the LDoS. 

As far as 7 is concerned the behavior of 2DW and the EBRM models are the 
same, and there is an indication of the crossover from the perturbative to the non- 
perturbative regime, as implied (in a non-rigorous fashion) by the LDoS theory. It 
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Fig. 17. The compensation time for 2DW simulations, and for the corresponding EBRM 
simulations as a function of T for various values of e. 
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Fig. 18. The estimated decay rate 7 for the same simulations as in the previous figure. 

is rather than 7 that exhibits sensitivity to the nature of the dynamics. This is be- 
cause is sensitive to the evolution of the main part of the distribution. We already 
had made this observation on the basis of the analysis of the WBRM model (see 
previous Subsection 8.1). Here we see another consequence of this observation. 



9 Conclusions 



There is a hierarchy of challenges in the study of quantum dynamics. The simple 
way to explain this hierarchy is as follows: Let us assume that there are two Hamil- 

tonians, Hi and 'H2, that differ slightly from each other. Let us then quantify the 
difference by a parameter e. Let us distinguish between a FOPT regime, Wigner 
regime, and non-perturbative (semicircle or semiclassical) regime according to the 
line shape of the LDoS. Do we have enough information to say something about 
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the dynamics? 



In the conventional wavepacket dynamics, one Hamiltonian is used for preparation 
and for measurement, while the other for propagation. It is well known that the 
Fourier transform of the LDoS gives the survival amplitude and hence V{t). But 
what about other features of the dynamics. What about the energy spreading 5E(t) 
for example? It turns out that the answer requires more than just knowing the LDoS. 
In particular we observe that in the non-perturbative regime physical models differ 
from the corresponding RMT model. In the former case we have ballistic spreading 
while in the latter we have diffusion. 

Is there any new ingredient in the study of driving reversal dynamics? Maybe it 
is just a variation on conventional wavepacket dynamics? The answer turns out to 
be interesting. There is a new ingredient in the analysis. This becomes very clear 
in the RMT analysis where we find a new time scale that distinguishes between 
a stage of "reversible diffusion" and a stage of "irreversible diffusion". This time 
scale (tsdn) can only be probed in a driving reversal experiment. It is absent in the 
study of conventional wavepacket dynamics. 

Things become more interesting, and even surprising, once we get into details. 
Let us summarize our main findings. We start with the conventional wavepacket 
dynamics, and then turn to the driving reversal scenario. 

The main observations regarding wavepacket dynamics are summarized by the di- 
agrams in Fig. 11. We always have an initial ballistic-like stage which is implied 
by FOPT. During this stage the first order (in-band) tails of the energy distribu- 
tion grow like t"^. We call this behavior "ballistic-like" because the second moment 
SE{t) grows like t"^. It is not a genuine ballistic behavior because the rth moment 
does not grow like but rather all the moments of this FOPT distribution grow like 
t\ 

The bandwidth is resolved at the time Td. In the perturbative regime this happens 
before the breakdown of perturbation theory, while in the non-perturbative regime 
the breakdown tpn happens before t^. As a result, in the non-perturbative regime 
we can get a non-trivial spreading behavior which turns out to be "ballistic" or 
"diffusive", depending on whether the system has a classical limit or is being RMT 
modeled. 

Once we consider a driving reversal scenario, it turns out to be important to mark 
the time tsdn when the energy distribution is resolved. The question is ill-defined 
in the perturbative regime because there the energy distribution is characterized 
by two energy scales (the "bandwidth" and the much smaller "core width"). But 
the question is well-defined in the non-perturbative regime where the distribution 
is characterized by one energy scale. It is not difficult to realize that for ballistic 
behavior t^dn ~ Td which is also the classical ergodic time. But for diffusion we get 
separation of time scales tp^ <^ ^sdn ^ 'Tch Thus we conclude that the diffusion has 
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two stages: One is reversible while the other is irreversible. 

But the second moment does not fully characterize the dynamics. In the other ex- 
treme we have the survival probability. Whereas SE{t) is dominated by the tails, 
V{t) is dominated by the "core" of the distribution. Therefore it becomes essential 
to distinguish between the FOPT regime where the "core" is just one level, and the 
rest of the perturbative regime (the "Wigner" regime) where the core is large (but 
still smaller compared with the bandwidth). 

The main findings regarding the driving reversal scenario are summarized by the 
following table: 



Regime 


perturbation strength 


PiT) behavior 


tr behavior 


SE{T) behavior 


1st order perturbative 


£ < £c 


LRT 

(super-Gaussian) 


LRT 


LRT 

(ballistic-like) 


extended perturbative 
("Wigner") 


£c < £ < £prt 


Wigner 

(Exponential) 


LRT(!) 


LRT 

(ballistic-Uke) 


non-perturbative 


£ > £prt 


non-perturbative 
(non-vmiversal) 


non-perturbative 
(non-vmiversal) 


non-perturbative* 
(diffusive/balUstic) 


*for the WBRM we have diffusion while for the 2DW model we have ballistic behaviour as impUed by classical LRT 



As expected we find that V(T) obeys FOPT behavior in the FOPT regime, which 
turns out to be super-Gaussian decay. In the Wigner regime 5E{T) still obeys LRT 
because the tails obey FOPT, while the non-perturbative core barely affects the 
second moment. But in contrast to that V{T) is sensitive to the core, and therefore 
we find Wigner (exponential) decay rather than FOPT (super-Gaussian) behavior. 
However, when we look more carefully at the whole V{t) curve, we find that this 
is not the whole story. We can characterize V{t) by the compensation time t,.. It 
turns out that is sensitive to the nature of the dynamics. Consequently it obeys 
"LRT scaling" rather than "Wigner scaling". This has further consequences that 
are related to quantal-classical correspondence. Just by looking at P(T) we cannot 
tell whether we look at the "real simulation" or on its RMT modeling. But looking 
on tr we can find a difference. It turns out that in the physical model tr exhibits e 
dependence, while in the case of RMT modeling tr is independent of e and exhibits 
"Wigner scaling". 

Finally we come to the non-perturbative regime. Here we have, in a sense a simpler 
situation. We have only one energy scale, and hence only one time scale, and there- 
fore 6E(T) and V{t) essentially obeys the same scaling. Indeed we have verified 
that the non-perturbative scaling with t^dn in WBRM simulations is valid for both 
the second moment and the survival probability. 
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Finally we would like to emphasize that the notion of "non-perturbative" behavior 
should not be confused with "non-linear" response. In case of quantized models, 
linear response of the energy spreading SE{T) is in fact a consequence of non- 
perturbative behavior. This should be contrasted with the WBRM model, where 
QCC does not apply, and indeed deviations from the linear response appear once 
we enter the non-perturbative regime. 

The study of irreversibility in a simple driving reversal scenario is an important step 
towards the understanding of irreversibility and dissipation in general. The analysis 
of dissipation reduces to the study of energy spreading for time dependent Hamil- 
tonians 7i{Q^ P]x{t)). In generic circumstances the rate of energy absorption is 
determined by a diffusion-dissipation relation: The long time process of dissipation 
is determined by the short time diffusion process. The latter is related to the fluctu- 
ations C{uj) via what we call "LRT formula". Thus the understanding of short time 
dynamics is the crucial step in establishing the validity of the fluctuation-dissipation 
relation. 
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